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ABSTRACT 

We present results from optical narrow- band (Ac = 8150 A and AA = 120 A) 
observations of the Great Observatories Origins Deep Survey (GOODS) fields, 
using Suprime-Cam on the Subaru Telescope. Using these narrow-band data, 
we then perform a survey of Lyman a Emitters (LAEs) at z ~ 5.7. The LAE 
survey covers an area of ~ 320 arcmin^ and a co-moving volume of ~ 8.0 x 10^ 
Mpc^. We found a total of 10 (GOODS-N) and 4 (GOODS-S) LAE candidates 
at 2; ~ 5.7. We perform a study of the spatial distribution, space density, and 
star formation properties of the LAEs at 2; ~ 5.7. 

Subject headings: cosmology: observations — cosmology: early universe galax- 
ies: formation — galaxies: evolution 



1. INTRODUCTION 

With the installation of the Advanced Camera for Surveys (ACS), the efficiency of the 
Hubble Space Telescope (HST) in performing wide-area and deep surveys has increased ten- 
fold. Such data are required to study formation and evolution of galaxies and are essential in 



^ Based in part on data collected at Subaru Telescope and obtained from the SMOKA science archive 
at Astronomical Data Analysis Center, which are operated by the National Astronomical Observatory of 
Japan. 

^Astronomical Institute, Graduate School of Science, Tohoku University, Aramaki, Aoba, Sendai 980- 
8578, Japan 

■^Space Telescope Science Institute 3700 San Martin Drive Baltimore MD 21218 
''^Also affiliated with the Space Sciences Department of the European Space Agency 
^National Astronomical Observatory, Mitaka, Tokyo 181-8588, Japan 
^INAF — Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy 



- 2 - 



identifying high-redshift galaxies. Significant progress has been made in recent years in per- 
forming multi- waveband deep surveys, greatfy extending tlie range of accessible redshifts and 
enhancing our knowledge of properties of galaxies out to the re-ionization epoch (Beckwith 
et al. 2004; Giavahsco et al. 2004a; Rix et al. 2004; Scoville et al. 2006 in preparation). 

One of the general methods to select galaxies at high redshift is drop-out technique (e.g., 

Steidel et al. 1996), monitoring the Lyman break (912 A) feature in galaxies as it moves 
with redshift to longer wavelengths and progressively out of the observed (shorter) passbands. 
Such high-z galaxies have been found in large numbers in the Great Observatories Origins 
Deep Survey fields (GOODS; Giavalisco ct al. 2004b; Dickinson et al. 2004; Stanway ct al. 
2004; Bouwens et al. 2004) and are used to investigate the origin of cosmic re-ionization at 
z > 6 (e.g., Barkana & Loeb 2001). Another technique to identify high- 2; objects is narrow- 
band search for strong Lyman a emitters (LAEs) (e.g., Hu, Cowie & McMahon 1998; Cowie 
& Hu 1998; Rhoads & Malhotra 2001; Hu et al. 2002; Ajiki et al. 2003; Kodaira et al. 
2003, Taniguchi et al. 2005). Such narrow-band surveys, conducted at selected wavelengths 
to avoid atmospheric emission, target emission line galaxies at certain redshifts. Combined 
with broad-band data to subtract the continuum, we can estimate the narrow-band flux for 
the respective lines observed. 

In this paper, we present results from a narrow-band survey in the GOODS fields, 
searching for LAEs at z = 5.7. The area coverage and depth of the GOODS, combined with 
its multi- wavelength imaging data, allows unbiased selection of a population of LAEs at high 
redshifts. Furthermore, availability of deep Spitzcr and Chandra data, high resolution ACS 
images and ground-based optical-IR data in the GOODS fields allow a study of the nature 
of LAEs, their morphology and star formation rate at z — 5.7 and provides an unbiased 
sample of LAEs to compare with homogeneously selected LBGs at a similar redshift. We 
adopt a flat universe with Jlmatter = 0.3, Qa — 0.7, and Hq —70 km s~^ Mpc"^. Throughout 
this paper, magnitudes are given in the AB system. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Narrow- band Observations and Source Detection 

Deep narrow-band surveys of the GOODS fields are performed using the Suprime-Cam 
(Miyazaki et al. 2002) on the 8.2 m Subaru Telescope. We use a narrowband filter, NB816, 
centered on 8150 A with a width AAfwhm = 120 A( Ajiki et al. 2003). This filter samples 
Lyman a emission lines in the redshift range 2; = 5.65 — 5.75. 

The NB816 observations of GOODS-S were made during 20-22 February 2004 (UT) 
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while, for the GOODS-N field, we used Subaru arcfiival data wfiicfi were taken in April 
2002 and 2003. All the NB816 observations were done under photometric conditions with a 
seeing between 0'.'7 and 0'.'9. The narrowband imaging data were reduced with the software 
package SDFRED (Yagi et al. 2002; Ouchi et al. 2004). The FWHM of stellar objects in 
the final NB816 images of the two fields are 0'.'9. The limiting NB816 magnitude in both 
fields is NB816 ~ 25.4 for a So" detection over a 2" diameter aperture. Deep broad-band 
images taken by the ACS (-B435, Vgoe; ^775) and 2:850) and ground-based near-IR data, taken 
as a part of the GOODS (Giavahsco et al. 2004a) are also available for these fields. 

The NB816 survey covers a total of 160 arcmin^ area in each field, corresponding to 
the full GOODS-ACS area (a total of 320 arcmin^ area). The transverse co-moving area of 
the LAE survey at ^ = 5.7 is 1.8 xlO^ Mpc^. Combined with the FWHM of the filter, which 
has a Gaussian-hke shape, this corresponds to a co- moving depth of 45 Mpc (zmin ~ 5.65 to 
-^max ~ 5.75) and a combine volume of 8.0 x 10^ Mpc^. 

Source detection was performed using SExtractor version 2.3 (Bertin & Arnouts 1996). 
A source is defined as 9-pixel contiguous area above the 3 a noise level on the NB816 image. 
Photometry was performed using IRAF - APPHOT with a 2" diameter aperture. To the 
magnitude limit of our survey, NB816 — 25.4 (f« 5(7 in the both fields), a total of 9,600 and 
8,700 sources are found in the GOODS-N and -S fields respectively. 



2.2. Selection of LAE Candidates 

The LAE candidates z ^ 5.7 are identified by combining the GOODS ACS images 
(S435, T/goe, ^775, and z^^q) (Giavalisco et al. 2004a) and narrow-band catalog with NB816 < 
25.4. The LAEs are selected to satisfy the following criteria: 

S435 > 27.5, (1) 
1^606 > 27.6, and (2) 
iz-NB816 > 0.7 (3) 

where iz is the continuum magnitude at A = 8150 A , estimated by linear interpolation 
between 2775 and ^850 flux densities {fiz = 0.63fi^^^^ +0.37fzgr,„)- The first two criteria ensure 
that the objects at z ^ 5.7 are undetected (at 2 a noise level) in 5435 and Vgoe while, the 
third condition allows selection of LAEs with a rest-frame Lyman a equivalent width > 17 
A. This is illustrated in Figure 1 where the LAE candidates are identified on iz — NB816 vs. 
NB816. color-magnitude diagram. A total of 10 (GOODS-N) and 4 (GOODS-S) LAEs at a 
targeted redshift of 2; ~ 5.7 are found to satisfy the above criteria. Their NB816 and broad- 
band ACS images are shown in Figure 2 (GOODS-N) and 3 (GOODS-S). Their coordinates 
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and photometric properties are summarized in Table 1. 

It is important to estimate the fraction of contaminants in our LAE survey. Prom 
spectroscopic observations of a sample of 23 LAE candidates at z ~ 5.7, identified in a 
narrow-band survey, Hu et al. (2004) found that 4 are interlopers [two are [O ll] and Ha 
emitters, one is a red star {NB816 — z' > 0.2) and one a source with no discernible spectral 
features]^. Since none of our LAE candidates have red NB816 — z^^q colors, similar to what 
expected for stars, our LAE survey is expected to be free from contamination by red stars. 
This, combined with the fact that the B- and y-band data here are deeper than those of Hu 
et al. (2004) by ~ 1 magnitude, we expect a contamination rate of <14% for our sample of 
LAEs here. 



3. Results 

3.1. Number Density and Lyman a Luminosities of LAE Candidates 

Using the estimated survey volume for each field, 4.0 x 10^ Mpc^, and the number 
of LAE candidates found at z 5.7, we estimate the LAE space density of n(LAE) ~ 
2.5 X lO-^Mpc-3 and ~ 1.0 x lO-^Mpc'^ for GOODS-N and GOODS-S respectively. The 
total space density in the combined GOODS fields is then n(LAE) ~ 1.7 x 10~^Mpc~^. The 
space density of LAEs in the GOODS-S is lower by a factor of 2.5 than that in GOODS-N. 
Spatial distribution of the LAE candidates are shown in Figure 4. 

In Table 2, we list the Lyman a luminosities for our sample of LAE candidates, as- 
suming them to be at z = 5.70, with the Lyman a shifted to the central wavelength of the 
NB816 filter. The Lyman a luminosities range from fti 4.2 x 10^^ ergs s~^ to 1.2 x 10^^ 
ergs s~^. These values correspond to - of the Lya luminosity function by Ajiki et al. 
(2003; = 4.1 ± 0.4 x 10^^ ergs s~^). In Figure 5, we show the distribution of Lyman a 
luminosities for our 14 LAE candidates (upper panel). We compare our results with those of 
the previous LAE surveys at z = 5.7 (Rhoads & Malhotra 2001; Ajiki et al. 2003; Hu et al. 
2004) (lower panel in Figure 5). It is clear that the number density of our sample of LAEs 
with Lyman a luminosities of 5 x 10^^ ergs s~^, is higher (by a factor of ~ 5) than those 
of the other surveys. This seems attributed to the relatively deeper completeness hmit for 
our survey [Liiin(Lya)~ 4.8 x 10^^ ergs s~^ corresponding to NB816 = 25.4] compared to the 
previous studies, including; Rhoads & Malhotra [2001; Liiin(LyQ;)~ 5.3 x 10*^^ ergs s~^], Ajiki 



^Hu et al. (2004) used the following criteria for selecting their LAE candidates; (a) NB816 < 25.05, (b) 
Ic - NB816 > 0.7, (c) undetected in the B and V images, and (d) r' - z' > 1.8. 
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et al. [2003; Lii^(Lya)c^ 7.0 x 10^^ ergs s-^], and Hu et al. [2004; Liii„(LyQ;)~ 6.6 x 10^^ 
ergs s~^]. The extra depth enables our survey to include more LAEs with the Lyman a 

luminosities ~ 4 6 x 10'*^ ergs s~^ than the previous surveys. We also note that, for 

luminosity [L(LyQ;)> 7 x 10^^ ergs s~^], the number density of our combined (GOODS-N 
and -S) sample is higher than that of Rhoads & Malhotra (2001) by a factor of «i 2 while 
lower than those of Ajiki et al.(2003) and Hu et al. (2004) by a similar factor. 



3.2. StEir Formation Rates 

Using the LAE luminosities, we now estimate their corresponding star formation rates 
(SFRs), using the relation 

SFR{Lya) = 9.1 x 10-^^L(LyQ;) Mq yr-\ (4) 

(Kennicutt 1998; Brocklehurst 1971). Here L(LyQ;) is in units of ergs s~^ and we assume a 
Salpeter initial mass function with (miowcD ""^uppcr) = (0.1 Mq, 100 Mq). The results are given 
in the third column of Table 2. We note that the SFR derived here can be underestimated 
due to the effect of absorption by H I gas both in the host galaxies and in the intergalactic 
space. The SFRs range from 4 to 11 Mq yr~^ with a median of 5.2 Mq yr~^. The estimated 
SFRs here are comparable to those of LAEs at 2; ~ 5.7 - 6.6 (e.g., Ajiki et al. 2003; Taniguchi 
et al. 2005). 

It is instructive to examine whether the SFRs derived from Lyman a luminosities here 
are consistent with those derived from the UV continuum luminosity from the broad-band 
data in the GOODS. We convert the observed ^350 magnitudes to UV continuum luminosities 
at A = 1360 A. Using the relation (Kennicutt 1998; see also Madau et al. 1998), 

SFR{\JV) = 1.4 X 10~2^L^ Mq yr^\ (5) 

where is in units of ergs Hz~^. The estimated SFRs are summarized in Table 2. Gom- 
parison between the SFRs for individual sources, estimated from SFR{Lya) and S'F_R(UV) 
(see Figure 6), shows that, on average, SFR(\JV) is relatively higher than SFR{Lya) for 
most of the LAE candidates. We also measure the average SFR ratios, iS'F/?totai(Lya) / 
S'Fi?totai(UV)= 0.71, where S'Fi?totai(LyQ;) and 5'Fi?totai(UV) arc, respectively, the sum of 
SFR{Lya) and S'Fi?(UV) of all our LAE candidates. This ratio is close to that obtained by 
Ajiki et al. (2003) for their sample of LAEs at z «i 5.7. The lower value of SFR{Lya) can 
be attributed to the effect of the dust extinction and scattering by the intergalactic medium. 
However, some of our LAE candidates have SFR{Lya) / SFR([JY) > 1. These are likely 
to be in a very early phase (< 10^ yr) of star formation activity, in which, SFR([JY) values 
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are underestimated (Schaerer 2000; see also Nagao et al. 2004, 2005). As illustrated in the 
right panel of Figure 6, objects with small SFRiXJV) have higher SFR(Lja) / SFR(\JY) 
ratios. This, in part, is explained by selection effects in the A'^fi^i^^limited sample (i.e., our 
survey cannot find objects below the dotted curve in the right panel of Figure 6). However, 
all objects with SFRiLya) / SFR{VY) > 1 (i.e., the age of < fO« yr) have SFR(UV) < 5. 
This implies that the UV continuum emission of very young star-forming galaxies are too 
faint to be found in typical LBG surveys. The narrow-band survey provides a very efficient 
method to investigate such very young star-forming galaxies. 

4. Summciry 

Wc present results from a NB816 survey in the two GOODS fields, aiming for LAE 
galaxies in the redshift range z = 5.65 — 5.75. We find a total of 14 LAE candidates satisfying 
our selection criteria. The space density of LAEs with L(LyQ;) 5 x 10^^ ergs s"-*^, found here 
(1.0 X 10~^ Mpc~^), is a factor of ~ 5 higher than other similar surveys (Rhoads & Malhotra 
2001; Ajiki et al. 2003; Hu et al. 2004), due to relatively deeper flux hmit in the present 
survey. Using the Lya luminosities, we estimate the SFRs for LAEs, SFRlLya) to be in the 
range 4 — 11 Mq yr~^. Overall, this is smaller than the SFR estimated from rest-frame UV 
light [SFR{Lja)/ SFR{Uy)= 0.71], due to extinction in Lya flux. We interpret a larger 
SFR{Lja) value [i.e. SFR(Lja) / SFR{l]'V)> 1], observed for a few LAEs as due to a 
very early phase of star formation (< 10* yr). 
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Table 1. the LAE candidates at ^ 5.7 in GOODS-N and GOODS-S fields. ^ 



No. 


a(J2000) 


5(J2000) 




Veoe^ 


1775^ 


2850*' 


NBSlff" 




h 


m s 


O 


f // 












GOODS-N 


1 


12 


35 46.4 


-1-62 


12 19 


29.26 


99 


26.91 


26.91 


25.12 


2 


12 


35 58.9 


-1-62 


10 18 


27.27 


27.87 


26.65 


26.81 


25.37 


3 


12 


36 07.9 


-1-62 


08 39 


27.56 


28.48 


26.17 


28.09 


24.79 


4 


12 


36 13.4 


-1-62 


07 48 


29.51 


34.20 


26.30 


25.44 


24.75 


5 


12 


36 45.4 


+62 


18 03 


29.07 


99 


28.35 


26.92 


25.40 


6 


12 


36 51.6 


-1-62 


19 37 


99 


28.53 


26.82 


25.76 


24.56 


7 


12 


37 01.0 


-1-62 


21 15 


99 


29.13 


26.10 


25.54 


25.03 


8 


12 


37 17.9 


+62 


17 59 


27.61 


31.58 


26.74 


26.17 


24.97 


9 


12 


37 27.8 


+62 


11 47 


99 


99 


27.43 


99 


25.23 


10 


12 


37 35.9 


+62 


14 22 


27.27 


29.08 


26.22 


25.45 


25.01 


GOODS-S 


11 


03 32 22.3 


-27 56 59 


27.34 


27.92 


28.40 


26.16 


25.26 


12 


03 32 29.4 


-27 40 29 


27.65 


27.24 


25.89 


26.48 


25.05 




03 32 37.. "5 


-27 40 57 


30.36 


28.93 


26.90 


26.87 


24.. 39 


11 


1):-! 


:-!2 52.1) 


-27 




28.1 J 


2!).()() 


2(i.<)!) 


!)<) 


2.5,17 



''AB magnitude in a 2" diameter. An entry of "99" indicates that no excess 
flux was measured 



Table 2. Lyman a luminosities and star formation rates for the LAE candidates at 2; ~ 5.7 



No. L(Ly«)=» SFR{Lya)'° Lvmo" SFR(UVf SFR{Lya)/SFR{VV) 

(10*2 ergs s-l) {Mq yr"!) (lO^S ergs s"! Hz"!) (Mq yr"!) 



1 


5.6 


5.0 


3.3 


4.6 


1.09 


2 


4.3 


3.9 


3.6 


5.1 


0.76 


3 


8.2 


7.4 


<2.0 


<2.8 


>2.65 


4 


6.5 


5.8 


12.7 


17.8 


0.33 


5 


4.2 


3.8 


3.3 


4.6 


0.83 


6 


8.7 


7.9 


9.6 


13.4 


0.59 


7 


4.7 


4.2 


11.7 


16.4 


0.26 


8 


6.0 


5.4 


6.5 


9.2 


0.59 


9 


5.6 


5.1 


<2.0 


<2.8 


>1.81 


10 


4.6 


4.2 


12.6 


17.7 


0.24 


11 


4.3 


3.9 


6.6 


9.2 


0.42 


12 


5.8 


5.2 


4.9 


6.9 


0.75 


13 


11.6 


10.5 


3.4 


4.8 


2.17 


14 


6.0 


5.4 


<2.0 


<2.8 


>1.92 



'^a Rj 0.1 X W*^ ergs s'^. 
t'o- « 1 M0 yr-i. 

^The UV continuum luminosity at A = 1360 A. a ^ 2.0 x 10^* ergs s'^ Hz"!. 
•^(T « 2.8M0 yr-l. 
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Fig. 1.— iz - NB816 vs. NB816 diagram in GOODS-N (left) and GOODS-S (right). Our 
LAE candidates are shown as filled circles. The vertical fines show the detection fimit corre- 
sponding to NB816— 25.4. The horizontal lines show the selection criterion iz—NB816 — 0.7. 
The curves show the 3(j limits for iz — NB816. 
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Fig. 2 — Broad-band and NB816 images of our LAE candidates at z 5.7 in GOODS-N. 
Each box is 12" on a side (north is up and east is left). Each circle has 3" radius. 




2. — continued. 
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Fig. 3. — Broad-band and NB816 images of our LAE candidates at 2; ~ 5.7 GOODS-S. 
Each box is 12" on a side (north is up and east is left). Each circle has 3" radius. 
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Fig. 4. — Spatial distributions of our LAE candidates in GOODS-N (left) and -S (right). 
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Fig. 5. — Upper panel: Distribution of Lyman a luminosity of our LAE candidates. Lower 
panel: Our results is compared with those from Rhoads & Malhotra (2001), Ajiki et al. 
(2003) and Hu et al. (2004). 
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Fig. 6. — Comparison between SFR(Lja) and SFR{l]V) (left panel) and SFR{Lya) / 
SFR{UV) and SFR{l]y) (right panel) for 14 LAE candidates. The dashed hnes in the 
both panels show SFR{Lya)/SFR{\jy) ~ 0.71. The dotted curve in the right panel show 
detection limit of our survey corresponding to SFR{Lya) — 3.8Mq yr~^. 



